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We study the out-of-equilibrium nonlinear dynamics of fields after post-inflationary preheating. 
During preheating, the energy in the homogeneous inflaton is exponentially rapidly transfered into 
highly occupied out-of-equilibrium inhomogeneous modes, which subsequently evolve towards equi- 
librium. The infrared modes excited during preheating evolve towards a saturated distribution long 
before thermalization completes. We compute the equation of state during and immediately after 
preheating. It rapidly evolves towards radiation domination long before the actual thermal equilib- 
rium is established. The exact time of this transition is a non-monotonic function of the coupling 
between the inflaton and the decay products, and it varies only very weakly (around 10~^^ s) as 
this coupling changes over several orders of magnitude. This result is applied to refine the relation 
between the number of efoldings N and the physical wavelength of perturbations generated during 
infiation. We also discuss the implications for the theory of modulated perturbations from preheat- 
ing. We finally argue that many questions of the thermal history of the universe should be addressed 
in terms of pre-thermalization, illustrating this point with a calculation of perturbative production 
of gravitinos immediately after chaotic inflation. We also highlight the effects of three-legs inflaton 
interactions on the dynamics of preheating and thermalization in an expanding universe. 



I. INTRODUCTION: BETWEEN INFLATION 
AND THERMALIZATION 

According to the inflationary scenario, the universe at 
early times expands quasi-exponentially in a vacuum-like 
state without entropy or particles. During this stage 
of inflation, all energy is contained in a classical slowly 
moving inflaton field (fi. The fundamental Lagrangian 
C{4>,Xt'4't ^ik, ■■■) contains the inflaton part with the 
potential V{(f>) and other fields which give subdominant 
contributions to gravity. In chaotic infiationary models, 
soon after the end of infiation the almost homogeneous 
infiaton field (j){t) coherently oscillates with a very large 
amplitude of the order of the Planck mass around the 
minimum of its potential. Due to the interactions of other 
fields with the infiaton in £, the infiaton field decays and 
transfers all of its energy to relativistic particles. If the 
creation of particles is sufficiently slow (for instance, if 
the infiaton is coupled only gravitationally to the matter 
fields) the decay products simultaneously interact with 
each other and come to a state of thermal equilibrium 
at the reheating temperature Tr- This gradual reheating 
can be treated with the perturbative theory of particle 
creation and thermalization as long as the couplings are 
sufficiently small Pj^jP]. However, the particle produc- 
tion from a coherently oscillating inflaton for a wide range 
of couplings occurs in a non-perturbative regime of para- 
metric excitation ^IHI^' This picture, with variation 
in its details, is extended to other infiationary models. 
For instance, in hybrid inflation inflaton decay proceeds 
via a tachyonic instability of the inhomogeneous modes 
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which accompany the symmetry breaking |3| . 

One consistent feature of preheating - non- 
perturbative copius particle production immediately 
after inflation - is that the process occurs far away from 
thermal equilibrium. The energy of the inflaton zero 
mode is transferred to particles in an out-of-equilibrium 
state with very large occupation numbers within a 
very short time interval of about 10"'^^ sec, which can 
be much shorter than the time needed for relaxation 
towards thermal equilibrium. 

One can expect the initial conditions for many parame- 
ters of the subsequent cosmological thermal history to be 
settled during or around preheating. It is often thought 
that the details of the transition between infiation and 
the hot radiation dominant stage are not relevant, except 
for the so-called reheating temperature T^, the tempera- 
ture of the ultra-relativistic plasma at the time when it 
reaches thermal equilibrium. Definitely, this is an impor- 
tant parameter of the early universe. However, a precise 
understanding of how thermal equilibrium is reached is 
crucial since partial thermal distributions can be respon- 
sible for cosmological baryo/leptogenesis, the possible 
creation of dangreous cosmological relics, etc. The out- 
of-equilibrium character of preheating opens the possibil- 
ity for equally relevant phenomena associated with non- 
equilibrium physics, including phase transitions, non- 
thermal production of heavy particles, etc. 

Precise knowledge of the expansion of the universe a (t) 
is also required to connect a physical wavelength of the 
observed cosmological perturbation, k/a (t), to the num- 
ber of e-foldings N at which that wavelength exited the 
horizon during infiation. 

Finally, we can have an alternative mechanism of gen- 
eration of (almost) scale free adiabatic metric perturba- 
tions from preheating, based on spatial modulation of 
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couplings Small spatial fluctuations in the cou- 

plings lead to small fluctuations of the rates of physical 
processes, which in turn generate adiabatic perturbations 
of the energy density and of the metric. Analyzing the 
details of this mechanism require precise knowledge of 
the evolution of the equation of state (EOS). 

To understand the early post-inflationary period we 
have to understand the dynamics of interacting fields in 
an out-of-cquilibrium state with large occupation num- 
bers evolving towards a state of ultimate thermal equilib- 
rium. This is a complicated problem of non-equilibrium 
quantum field theory, which is by itself a very interesting 
topic. 

The theory of the transition from inflation to thermal- 
ization has been investigated with fully non-linear nu- 
merical lattice simulations plus different techniques of 
classical field dynamics jllllll [H Classical field 

dynamics is adequate as long as occupation numbers of 
field excitations are large. There has also been progress 
in understanding out-of-equilibrium QFT dynamics in 
0(N) sigma models beyond the Hartree approximation 
|l5|. taking into account crucial effects of rescattering 
after preheating From all of these methods, the fol- 
lowing picture emerges. Immediately after preheating, ei- 
ther from parametric resonance or tachyonic, one or more 
bose fields are excited in an out-of-equlibrium state with 
large occupation number. Backreaction of these fields 
terminates their production. Interaction (rescattering) 
of these modes between each other, and with the remain- 
ing infiaton field after the first stage of preheating is vi- 
olent and non-perturbative. During this very short stage 
a large amount of entropy is generated, chaotic (turbu- 
lent) wave dynamics is established, and the fields not 
directly excited during preheating get excited in out-of- 
equilibrium states due to rescattering; residual infiaton 
oscillations are still present. A next, longer stage then 
takes place, characterized by lower (but still large) oc- 
cupation numbers such that rescattering is perturbative, 
and the occupation spectra gradually move towards a sat- 
urated state by cascading towards ultra-violet (UV) and 
infra-red (IR) modes (in the spirit of Kolmogorov wave 
turbulence). The last and the longest stage will be the 
stage of proper thermalization, when the distributions 
evolve towards thermal equilibrium. Quantum physics 
is important at the end of this stage, when the occupa- 
tion numbers are small and the classical approximation 
breaks down. One may say that reheating is completed 
when this last, longest stage is completed, and its timing 
defines Tr and a^. 

In this paper we suggest a new look at the transition 
between infiation and thermalization. We focus on cal- 
culation and understanding of the effective equation of 
state w = p/e throughout all of these stages after infia- 
tion. We use numerical lattice simulations to cal- 
culate w in the simple chaotic model with a quadratic 
potential and a simple four-legs interaction g^cf^y^ . 

In particular, we notice that w approaches (but 
does not necessarily reach) the equilibrium radiation- 



dominated value 1/3 while the fields are far from ulti- 
mate thermal equilibrium. We compare this with the 
analysis of devoted to thermalization of the O {N) 
sigma model in Minkowski spacetime. In that case, the 
equation of state was found to evolve sharply towards 
w — \/Z long before thermalization completes, which 
prompted the authors of [T^ to describe the state of the 
system as a pre-thermalized state. Although we also see 
a somewhat similar effect, the expansion of the universe 
and the presence of a residual massive infiaton compo- 
nent, prevent the equation of state from being exactly 
the one of radiation at this very early stage. Moreover, 
we enphasize that the exact microphysics of the system 
is described by a turbulent state, which is still very far 
from a thermalized one. (This difference can be impor- 
tant when discussing the cosmological effects of decay 
products, as for instance production of dangerous relics 
from nonthermal distributions.) Indeed, our simulations 
show that the occupation numbers of excited infra-red 
modes evolve towards a steady state, related to turbu- 
lence of classical interacting waves. 

The plan of the paper is the following: 

In Section^we describe preheating in a simple chaotic 
infiationary model. In Section IlIII we present the results 
of lattice simulations for time-dependent variables, in- 
cluding the equation of state w(i) during different stages 
of preheating, occupation numbers and fiuctuations. 

In Section im we discuss some immediate applications 
of our results for cosmology. We consider the application 
of the EOS evolution to modulated cosmological fiuctua- 
tions and to the N — log k formula. We note several qual- 
itatively important issues in the dynamics of thermaliza- 
tion. The interaction (jy^ describes (at the perturbative 
level) scattering between infiaton quanta, rather than in- 
fiaton decay. Scattering between massive </> soon becomes 
inefficient, and does not lead to a complete depletion of 
the infiaton, which eventually ends up dominating over 
the light degrees of freedom x- This problem is automat- 
ically avoided if three- legs bosonic interactions (^ly^ are 
also present. We illustrate this with an example moti- 
vated by supersymmetry. We also discuss the rapid sat- 
uration of the IR modes excited during preheating, which 
occurs long before the ultimate thermal equilibrium (of 
all modes) is reached. 

One consequence of the excitation of infrared modes is 
that it can lead to an overproduction of gravitinos or of 
other dangerous gravitational relics. This generation is 
different from direct nonperturbative production during 
preheating Instead, it is analogous to the pertur- 

bative thermal production first considered in "IQ^. The 
main difference with [191] is that we compute the produc- 
tion right after the first preheating stag e, and not only 
after thermalization has completed MJ]. Although the 
distributions formed at preheating/rescattering are far 
from thermal, the high particle occupation numbers can 
lead to a significant generation of gravitinos. 

The above results are summarized in Section Ivl 
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II. THE MODELS AND METHODS 

In this section we briefly describe the model of chaotic 
inflation with a quadratic potential and additional inter- 
action terms and the methods we use to investigate the 
fields dynamics. The quadratic form of the potential is 
a good generic approximation around the minimum. For 
simplicity, we take the potential 

V{<P) = im202 

also throughout inflation. The inflaton mass is taken to 
be m ~ 10^^ Gev, which, in the simplest cases, repro- 
duces the correct amplitude for the CMB anisotropics. 
After inflation, the background inflaton field oscillates 
with its amplitude diluted by the expansion as </> oc 
(/>o/a^/^, where 0o ~ Mp/lO is the initial amplitude, 
and a{t) is the scale factor of the FRW flat universe. 
More accurately, the inflaton evolution is given by (f)(t) « 
—M^sinmt. The energy density of the inflaton oscil- 

lations evolves as e w m^c^g/a'^, while the pressure is 
p w — (m^0g/a^) cos2mi. Since the frequency of the 
oscillations is higher than the rate of expansion H, we 
deal with a nearly vanishing averaged effective pressure, 
p « 0. Therefore, immediately after inflation, the equa- 
tion of state (EOS) is effectively the one of matter dom- 
ination 

w=p/€^0. (2) 

This is what we expect for inflaton oscillations inter- 
preted as the coherent superposition of heavy inflaton 
(quasi)-particles at rest. 

Inflaton oscillations decay due to the coupling between 
the inflaton and other fields. One may think about the 
inflaton coupling to other bosons through four-legs or 
three-legs interactions, for instance coupling to another 
scalar x as (?^</>^X^ or Yukawa-type couplings to 

fermions will be of the form hip^'ip. During preheating 
the leading channel of inflaton decay will be parametric 
resonance decay into bosons. A three legs interaction 
between bosons g^uipx^ requires a dimensionful scale cr, 
which emerges naturally in SUSY theories. However, in 
most cases the amplitude of the inflaton is sufficiently 
high during the early stages of preheating that the four- 
legs interaction dominates |2^. For this reason, at this 
stage we consider only the four-legs interaction 

A„t = -y-^'x' • (3) 

Three-legs interactions are important for the completion 
of the decay of massive inflatons, and we will discuss 
them further in Section HVI and in detail in [20l |. 

Before proceeding with the generation of x particles 
from the inflaton oscillations, we need to comment on the 
range of the coupling . At the level of QFT we can ob- 
tain an upper limit for . Indeed, in non supersymmetric 



models, the interaction Q leads to radiative corrections ^ 

to the effective potential ((TJ , of the form 12^5' l'^ 4'- 
inflaton values (p ~ 4Mp , which correspond to the scales 
where cosmological fluctuations are observed, radiative 
corrections do not alter the the potential ^ as long as 
g^ ^ 10~^. In supersymmetric theories radiative correc- 
tions from bosons and fermions cancel each other out so 
that the value of g can be even higher. 

There is also a lower limit on g^. During inflation 
the field x has the effective mass -I- 5^(0^) ~ g^Mp , 
where is its "bare mass" : we take <^ m and thus 
neglect it. It is known 21] that if we have two scalar 
fields (p and Xj the latest stage of inflation will be driven 
by the lightest scalar (light scalars are always generated 
during inflation). A consistent setting where the end of 
inflation is driven by the (j) field requires g'^M^ > m^, i.e. 
g^ > 10~^^. This lower limit can be evaded, however, if 
X is conformally coupled to gravity. 

A quantum scalar held x in a flat FRW background 
has the eigenfunctions Xk{'t)e-~^^^ with comoving mo- 
menta k. The temporal part of the eigcnfunction Xfc(^) = 
a^^'^Xk{t) obeys the equation 

Xk+[^+g^<f^Xk^Q (4) 

with positive- frequency initial conditions in the far past. 
In the effective frequency ux^ « ^ + 3^'/'^ we neglect 
small terms proportional to H^. This is an oscillator-like 
equation with a periodic frequency, which contains expo- 
nentially unstable solutions corresponding to parametric 
resonance. It is convenient to use a new time variable 
z — mt. Then, the essential dimensionless coupling pa- 
rameter which characterizes parametric amplification is 

2,2 

q — . The limits on the coupling g we have discussed 
above correspond to 10^^ ^ q ^ 10^. A large portion of 
this interval is in the broad parametric resonance regime 
{q > 1) of production of x particles 0. 

What is important, however, is that even one bosonic 
field with a large parameter q is sufficient to have the 
leading channel of inflaton decay occur through broad 
parametric resonance. Indeed, if one such field x is ex- 
ponentially amplified in this way then other flelds that 
are coupled to it directly or indirectly will tend to be 
amplified as well [T^ . 

We are interested in the dynamics of the flelds cj) and 
X after inflation. Although this toy model is simple, its 
dynamics are rather rich and complex. The first stage is 
the stage of parametric excitation of the quantum fluc- 
tuations of the field x, where the light x particles are 
copiously produced from the classical inflaton oscillation. 



There is also a divergent quadratic contribution. However, for 
the supersymmetric theory this term is independent of the value 
of (f) and can thus be cancelled by a counter-term. 



4 



This process is described by the theory of broad stochas- 
tic parametric resonance in an expanding universe Q. 
As soon as the occupation number of created particles 

"'= = t1^^ + "^'=" )-2 

becomes large, Uk S> 1, one can approximate the dy- 
namics with that of the classical interacting fields </> and 
X- Fully non-linear wave equations for coupling clas- 
sical fields are rather complicated. We use the LAT- 
TICEASY numerical simulations 0| to solve the equa- 
tions Dcf) + m'^(j) + g^4'x'^ = and + g'^4>'^X = in 
an expanding universe (□ being the covariant wave op- 
erator in the FRW geometry). The scale factor a{t) is 
self-consistently calculated from the Friedmann equation. 
The initial conditions correspond to gaussian initial fluc- 
tuations of (j) and Xi which arise from vacuum fluctua- 
tions. 

The novel element of our current study is the focus on 
new diagnostics of the numerical experiments. In par- 
ticular, we investigate the evolution of the equation of 
state w{t) during different stages, and for various cou- 
plings . Many previous numerical studies of preheat- 
ing concentrated on inflation with a quartic potential 
1^(0) = \4>^ 1 where the expansion of the universe can 
be "scaled away" . The inflaton oscillations in this model 
have the EOS w = 1/3, which corresponds to that of 
radiation domination, so w does not change after the in- 
flaton decay. By contrast, w changes during preheating 
in the more general model with a quadratic potential. 

The total energy density e and pressure p contain con- 
tributions from kinetic and gradient energies of the fields 
(j) and Xi plus interaction terms. It is instructive to cal- 
culate separate contributions from all of these terms. We 
find that the interaction term, except for a transient mo- 
ment between the end of preheating and onset of turbu- 
lent regime, is sufficiently small that we can interpret the 
results in terms of different particles 54> and x coupled to 
each other. In particular, this allows us to interpret the 
energy density in each component in terms of the integral 

where the comoving frequency is 



t^fc = yfc^ + JTi^Z/O^ , (7) 

and the co- moving occupation number Uk is given in jSJ. 
To study w we have to monitor the relationship between 
the effective masses of the fields 

"^le// ="^'+5'(x') . (8) 



ml,ff=g^{<Pf+g\5^^). (9) 



and the typical momenta. 

The variances of the fluctuations 

are also an important diagnostic of the non-equilibrium 
physics. Let us rewrite the wave function Xk{t) in 
the WKB form Xkit) = g-*/ ''*'^'= + dtu^^ 
While the occupation numbers in terms of a — /3 coeffi- 
cients are simply Uk = |/3fcP, the exact expression for the 
variances is more complicated 

<-^'>^^/f + *•"'))■ 

(12) 

Particles x produced at preheating are created in a 
squeezed state. In particular, this is reflected in the pres- 
ence of the oscillatory second term in the integral l|12|l . 
where the phases of a^, (3k are present. 

The total co-moving number densities of x a-nd 4> par- 
ticles are 

^-=(2^/^'^"- (13) 

^^-(^/rf^^-t, (14) 

respectively. It is instructive to monitor the time evolu- 
tion of the total number of particles A'tot — + N^p. At 
leading order in g'^ the total number of particles is con- 
served, A^tot (t) = const. However, there is a brief period 
immediately after preheating where A'tot is poorly defined 
due to the strong interactions between the waves. Indeed, 
during preheating the occupation number of x particles 
can be as large as Uk ~ l/ff^- Consider the collision in- 
tegral in the kinetic equation for the particle. If we for- 
mally represent it as a series with respect to the coupling 
5^, the smallness of will be compensated by the large 
occupation numbers. This means that the next terms of 
the collision integral associated with the higher order dia- 
grams (for example, the diagram with four incoming, two 
outgoing legs, and two vertices) are as important as the 
first term. In fact, the non-perturbative character of the 
interactions during this short intermediate stage makes 
it difficult to investigate the whole dynamics analytically. 
After this period ends, the spectrum of the occupation 
numbers formed at preheating is weighted towards the 
infra-red relative to a thermal distribution of the same 
energy density. Hence, particle fusion is more favored 
than particle dissociation. Once rescattering leads to a 
decrease of n^, only the leading four- legs diagram with 
one vertex becomes important, and A^tot decreases on 
the way to thermalization, on a timescale which is much 
longer than the one which can be observed in our lattice 
simulations. 
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III. OUTPUT OF THE CALCULATIONS 
A. The Calculations 

We performed three-dimensional lattice simulations for 
the model of Section |n] Our grid was a 256 x 256 x 256 
cube with a comoving edge size L = lO/m, which corre- 
sponds to a comoving grid spacing of dx ~ 0.04/to. As 
energy flows towards the UV end of the spectrum the 
simulations eventually reach a point where the grid spac- 
ing is too large to capture the important UV physics. 
By monitoring the spectra of the fields, however, we can 
verify that these simulation parameters were adequate to 
capture the relevant IR and UV physics well past the end 
of preheating. The time step was dt = 0.001/m and the 
inflaton mass was m = 10"*^ Mp. We used values of the 
coupling near = 10~^ This value is optimal because 
it is large enough to produce highly efficient preheating, 
but small enough that the occupation numbers Uf^ 1 j 
produce strong rescattering. The results should be quali- 
tatively similar for a wide range of values of g^, but would 
require more IR and/or more UV to simulate accurately. 

To probe later times and wider ranges of the couplings 
it will be necessary to extend the lattice simulations. 
This can be done with a parallelized version of the sim- 
ulation code LATTICEEASY (currently under construc- 
tion), or by combining the straightforward lattice simu- 
lations with other methods, like the equations for a large 
number of weakly coupled oscillators [221 ■ We intend to 
pursue both of these approaches in subequent work. 

In the rest of this section wc present the results of our 
simulations. 



B. Equation of State 

The time evolution of the EOS w[t) for different cou- 
plings is shown in Figure ^ Each point plotted on this 
figure represents the value of w averaged over a complete 
inflaton oscillation. This represents one of the main re- 
sults of our study. 

Immediately after inflation, the EOS averaged over in- 
flaton oscillations is w = 0. It sharply changes at the end 
of preheating. 

There are at least three important points worth enpha- 
sizing about the evolution of w. 

i) First, the transition of the EOS from w = to the 
value w 0.2 — 0.3 occurs very sharply, within a time 
interval ~ 10^'^^ sec. 

Indeed, recall that the unit of time on the plots is 1 /m, 
where m is the inflaton mass, i.e. 10""^^ sec. The first 
stage of preheating is completed within about a hundred 
of these units, i.e., 10"'^^ sec. The rise of w and gradual 
saturation takes roughly the same time. 

ii) Second, the dependence of w{t) on the coupling g^ 
for resonant preheating is a non-monotonic function of 
9^- 
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FIG. 1: Evolution of the equation of state w = w{t) as a 
function of time (given in units of m~^) for various couplings 
g'^ around = 2 x 10"'^. 

This is to say that the time during which preheating 
comes to an end is very weakly (logarithmically) depen- 
dent on the coupling. As seen from Figure ^ the curves 
w{t) begin to shift to the left towards an earlier end of 
preheating, as we vary g^ by 5%. However, at some point 
the curves stop moving to the left and instead begin to 
return toward the right. As we change g^ by about 25%, 
the cycle repeats. As we vary g^, the function w not 
only shifts, but it also varies its detailed shape. Still, 
to characterize these variations, we pick up the moment 
where w is equal to the value 0.15 (just for convenience 
of calculation), w{ttran) — 0.15. This allows us to plot 
the transition moment ttranig^) s& a function of 5^, see 
Figure |21 
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FIG. 2: Transition (preheating) time as a function of . 

We see that the transition time varies between 100/m 
and 150/m. This non-monotonic behavior of the dura- 
tion of preheating is explained in the theory of broad 
paremetric resonance @ (see Sections 6 and 9 there). 

The g^ dependence of the EOS is the critical issue 
for the theory of modulated cosmological perturbations, 
which we will discuss in Section llVI 
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iii) The third point is that w does not necessarily im- 
mediately go to the radiation dominated value 1/3. This 
is partly because immediately after preheating the light 
field stiU has a significant induced effective mass due to 
the interaction, and partly due to the significant residual 
contribution from the homogeneous inflaton Unfor- 
tunately, hmitations on running longer simulations pre- 
clude us from seeing further details of the time evolution 
of w. However, we have a strong theoretical argument to 
advance the discussion further. In a model with a massive 
inflaton and light scalar x even the radiation dominated 
stage is transient. Indeed, sooner or later the massive 
inflaton particles, even if significantly under-abundant at 
the end of preheating, will become the dominant compo- 
nent, and the universe will again be matter-dominated. 



C. Occupation Numbers 

The occupation numbers and nf. are among the 
most interesting variables to understand the micro- 
physics in our system of two interacting fields. First, 
we shall determine when the definition of Uk is mean- 
ingful. To do so, we consider the composition of total 
energy density e of the system of coupled (p and x fields. 
The total energy density e can be decomposed into par- 
tial contributions from the kinetic energy of both fields, 
their gradient energy, their potential energy (0 only in 
this model), and finally the interaction energy 
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(15) 

Figure 131 shows the relative contribution to the total en- 
ergy from each of the components for — 2.5 x 10^^. 



© for the energy of the particles is not a good approxi- 
mation and the occupation number Uk is not well defined. 
Outside this time interval, Uk is a meaningful quantity. 
Secondly, the contribution from the background homoge- 
neous inflaton is dominant even after preheating, up to 
approximately t ~ 150/m. A similar point was made for 
the model in 0. 

Let us now turn to the occupation numbers Uk- We 
find it more instructive to output not the occupation 
numbers Uk per se (as it is commonly done in the lit- 
erature) but the combination of Uk with the energy per 
mode ujk ■ This combination can be immediately com- 
pared with the Rayleigh- Jeans spectrum, 

n, « , (16) 

UJk- fi 

which corresponds to the equipartition spectrum of clas- 
sical waves (we introduce the chemical potential for 
generality). The comparison allows to determine how 
close the distribution is to the thermal one. 
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FIG. 3; Relative contribution of each of the energy compo- 
nents to the total energy, as a function of time. The vertical 
axis is the log of the various energy components in units of 
the initial energy m^^g multiplied by a^. 

We note two features of this plot. First, the interaction 
term is comparable to the other terms in the time interval 
100/to < t < 120/m. In this short period, the formula 



FIG. 4: Time evolution of the combination Lo^^knjl, for the 
model = 2.5 • 10"^ 

The combination Uk i^k for the two fields is shown at 
some characteristic times in the two Figures 0] and 
There are three distinct stages which characterize the 
evolution covered by our simulation. The first stage, 
characterized by linear dynamics, is the one of pre- 
heating and early rescattering. The first modes to be 
populated are the IR ones. Preheating of x particles 
occurs in the resonant band at comoving momentum 
fc» = y/gm(f>Q a^^^ ~ 7ma^^'^ 6J. Then, quanta S cj} are 
generated by rescattering. The annihilation Sxk Sxk — *■ 
5(l)k S(f)k amplifies quanta of the inflaton at A: ~ fc* . Even 
more effective is the rescattering of the x quanta against 
the infiaton zero mode, Sxk 4>o ~^ ^Xk S(j)k , which pro- 
duces quanta of both and x &t momentum fc ~ A:*/2 . 

The second stage is a violent stage of highly nonlinear 
dynamics. Starting at mt ~ 100 — 110 , the higher band 
at A:*/2 both increases in its amplitude and broadens to- 
wards higher momenta; quite interestingly, the peak lo- 
cation shifts from ~ fc*/2 to ~ during this quick and 
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FIG. 5: Same as Fig.H but for the x field. 



explosive stage. As we remarked, the particle occupation 
number is ill defined at this stage, and a nonlinear wave 
description is more adequate. 

The explosive stage of rescattering ends at about 
mt ~ 130 . In the next stage, characterized by pertur- 
bative dynamics, the distributions smooth out and start 
evolving (at a much slower rate) towards higher comov- 
ing momenta. The spectra in the IR approach a satu- 
rated power-law state, which then slowly propagates to- 
wards the UV. Although one can observe a greater ten- 
dency towards thermalization for the IR modes (where 
the rescaled spectra are closer to flat), the overall dis- 
tributions are still typical of the turbulent regime, and 
they are far from thermal. If we contrast this with the 
macroscopic behavior we have described above, particu- 
larly the evolution of the EOS, we see that the system can 
be considered in a pre-thermalized state (see also \vi\). 
but that thermalization is still far from being complete . 

It is also instructive to consider the product of the 
occupation number with the phase space sphere area 

and energy per mode ujk ■ This combination represents 
the energy density of the quanta at momentum k (since 
the product 4:Tt k"^ uij-nj-dk is the energy density of the 
quanta whose momentum has a magnitude between k 
and k + dk). 

In Figures IHl and [7| we plots the distributions at dif- 
ferent times separated by intervals St = An/m. This 
allows us to visualize the growth of the distributions, 
and to monitor the cascade of energy in the phase space 
k. When the ultraviolet part of the distribution hits the 
highest momentum of the simulation (defined by the grid 
size), the energy is artificially refiected back to the IR 
modes and the simulation is no longer reliable. Due to 
the scale chosen (natural rather than log scale) , only the 
distributions at times greater than about 100/to can be 
appreciated in the two figures and |7| Moreover, the 
double peak structure that can be observed for x at late 
times is due to the rescaling chosen. The plot for the 
occupation number, Uk k^ (not shown) has a high peak 
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FIG. 6: Time evolution of the combination k uikUk for the 
field (j) . The thicker (green) curve is the spectrum at the final 
time of our evolution. 




FIG. 7: Same as Figure El but for the field X- 



at A: ^ 4m , followed by a plateau up to fc ~ 30 m . The 
combination k'^uj^n^ shown has a peak at the momenta 
corresponding to this plateau, showing that these mo- 
menta dominate the energy density in the x distribution. 
The saturated spectra can be clearly seen in the green 
(thick) curves on these plots. 

The combination k'^tUkUk is also given in Figures|Hland 
13 where we however show only few times, and we focus on 
the IR part of the distributions. These figures show the 
rapid broadening (towards the UV) of the distributions 
in the violent rescattering stage, and the tendency to 
saturation at later times. 



D. Fluctuations and Effective Masses 

The evolution of the scalar fields can be strongly af- 
fected by the presence of dynamical effective masses for 
the two fields. A high mass for some of the fields can 
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of the variances (5</>^) and (x^)- Fields squared are rescaled 
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FIG. 11: Evolution of the scale factor a{t). 

FIG. 9: Same as Figure |H1 but for the field X- 



disfavour or block some of the interactions, with a con- 
sequent delay of thermahzation. For instance, it is weU 
known that a condensate which is strongly coupled to 
some field "0 can have a very long lifetime, since the same 
coupling which would allow the decay also produces a 
high effective mass for -0, which can prevent the decay 
from occuring. In the case we are discussing, the effec- 
tive masses of the fields acquire significant (loop) con- 
tributions from the high variances produced at preheat- 
ing/rescattering, according to Eqs. JHl-®- 

Figure [TUl shows the evolution of the background infla- 
ton oscillations and fluctuations of (50^) and (Sx^)- All 
values are plotted in "comoving scales", i.e. multiplied 
by a^. The corresponding physical quantities can be ob- 
tained by the value of the scale factor, which is shown in 
Figure ^ 

Figure^] shows the effective masses for and x i mul- 
tiplied by the scale factor a. This combination enters in 
the comoving dispersion relations for the quanta of the 
two fields, given by eq. 0. Therefore, the product arrics 



is the correct quantity to compare with the comoving mo- 
mentum k . 

This comparison is shown in Figure 1131 where we 
show the fraction of relativistic quanta (that is, with 
k > amcs) for the two fields. We see that the majority 
of the particles are non-relativistic, so that the effective 
masses indeed play a significant role in the evolution of 
the spectra. In particular, we observe that during rescat- 
tering very few particles are relativistic; this confirms the 
fact that the distributions are very peaked in the IR at 
this stage. The fraction of relativistic x particles rapidly 
increases during the violent thermahzation stage. The 
fraction of relativistic (50 remains instead always smaller 
(confirming the fact that m^ off > ^x,cS )i s-nd if slowly 
increases during the thermahzation stage, when the dis- 
tributions evolve towards the UV. 

Up to mt ^ 120 , the coherent mode of the inflaton 
dominates the effective mass for x ■ (The rapid variation 
of m-^^eg due to the coherent oscillations of is the cause 
of preheating.) At later times, the variance of gives in- 
stead a greater contribution. In the effective mass for 
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FIG. 13: Fraction of particles in the relativistic regime, k > 

the "bare" term m dominates up to mt ~ 100 , while the 
variance of x provides the dominant contribution at later 
times. Eq. (|12|l provides a useful expression for the vari- 
ance; during the turbulence stage, the occupation num- 
bers vary only adiabatically, so one may expect that the 
growth of a rrice results in a decrease of the variance, and 
vice-versa. This relation is visible in the figures shown. 
The "late time" decrease of a m-^^cS is accompanied by 
an increase of (x^)- This in turns cause an increase of 
cLm^^cS, and, consequently, a decrease of {(j)'^)- 

IV. COSMOLOGICAL IMPLEMENTATIONS OF 
THE RESULTS 

A. Complete Decay of the Inflaton 

In the simple model discussed so far, with the inter- 
action term (j)^ i the decay of the inflaton is not 
complete. Indeed, such a term mediates scatterings be- 
tween inflaton quanta, rather than single particle decay. 



Once (f> is diluted by the expansion of the universe, the 
scatterings become inefficient, and the number of infla- 
ton quanta remain practically constant. This poses a 
problem because in order to reach the stage of radiation 
domination we need to have complete decay of the infla- 
ton. 

The simplest way to have a complete decay is to con- 
sider three legs interactions. For instance, we can re- 
place © by 

Cint = -^ {cP + afx' (17) 

where a is a mass dimension parameter, which breaks 
the (j) — (j) symmetry. For a <^ <j} , the trilinear term 
produced by this interaction is irrelevant and the infla- 
ton decay occurs as described above. However, when 
is decreased by the expansion of the universe, the trilin- 
ear interaction becomes dominant, leading to a complete 
inflaton decay. 

We can avoid discussing the introduction of the scale 
a by considering a Yukawa interaction of the inflaton 
with some fermions tp . By itself, this interaction has 
interesting consequences for preheating, as discussed e.g. 
in ■ When both fermionic and bosonic interactions are 
present, preheating into bosons is typically much more 
efficient, since the one into fermions is reduced by Pauli 
blocking. However, once the interactions become per- 
turbative, the trilinear (jj ipip decay will eventually 
dominate over the 4>4> — > XX scattering. 

Three-legs interactions arise very naturally in SUSY 
theories. Indeed, consider the simple superpotential 

W^=^02^gW- (18) 

The corresponding scalar field potential contains three- 
legs as well as four-legs plus self-interaction terms 

y = + V2g'm^x' + ^gW + g'x' m 

(neglecting the imaginary parts of the fields). The four- 
legs interaction dominates over the three-legs one as long 
as >, (T = 2g^m. For reasonable values of a, this hap- 
pens during the first stages of preheating. Eventually, the 
amplitude of 4> decreases due to the expansion of the uni- 
verse, and the trinlinear interaction dominates, resulting 
in a complete decay of the massive inflaton. 

Hence, considering a supersymmetric theory automat- 
ically introduces trilinear vertices among the scalars. It 
also allows for stronger couplings between the inflaton 
and other fields, without spoiling the flatness of the in- 
flaton potential, as we have discussed in Sec. 

B. Relation between efoldings A'^ and wavelength of 
perturbations 

The precise history of the expansion of the universe 
a (t) is important for connecting the physical wavelength 
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k/a{t) of the cosmological perturbations we observe to 
the number of efoldings N at which the perturbation was 
generated (left the horizon) during inflation. 
According to the general lore [24. l25j | , one has 



N{k) = 62 - In 



k 



6.96 X 10-5 Mpc" 



A, (20) 



where 6.96 x 10~5Mpc~^ is the inverse size of the present 
cosmological horizon and A is defined by the physics after 
inflation: 



A 



-In- 



IQi^GeV 



1/4 



1, Vk 
-m 

4 Vend 



1 v^'t 

_j _cnd_ 
3 1/4 
/^rch 



(21) 



where Preh is the energy density at the end of reheating, 
Vfe is the value of the inflaton potential at the moment 
when the mode with the comoving wave number fc exits 
the horizon at inflation and finally T4nd is the value at 
the end of inflation. 

The last term in the right hand side of eq. (|21|l is re- 
lated to reheating, specifically that it is not an instan- 
taneous process. It is typically computed under the as- 
sumption that the universe is matter dominated (by the 
coherent oscillations of the inflaton) all throughout re- 
heating, and that it (instantaneously) becomes radiation 
dominated at the end of the reheating stage. As we have 
seen, this is not what happens. 

Indeed, a better definition of the last term in A is 
— iln-^ where Ord is the scale factor at the moment 
after reheating at which the EOS reaches its value at 
radiation dominance. As we have been discussing, this 
moment can occur while the distributions of particles are 
very far from thermal, shortly after the end of preheat- 
ing. This amounts to taking ard/oinf ^ 10 — 20, rather 
than Ord/ainf ^ 10* , where ord is the value of the scale 
factor at the end of reheating under the common assump- 
tion that perturbative reheating occurs entirely within a 
regime of matter domination, with T^. ~ 10^ GeV. In this 

situation one has i?RD ^ at the end of matter dom- 
ination, and the corresponding value of the scale factor 

/ v2/3 / v2/3 

ORD can be estimated as - Vh^) ~ ( j • 
Hence, for a preheating scenario with a very wide range 
of values of the coupling , the first term of can be 
fixed to 



iln^^_(0.6-0. 

4 flinf 



(22) 



This result holds for a wide class of infiationary models, 
and for many different types of preheating. 



plasma after the infiaton field decays and transfer its en- 
ergy into radiation.^ This is traditionally called the re- 
heating temperature Tr. As we have seen in many pre- 
vious papers and above in this paper, ultimate thermal 
equilibrium is preceeded by several stages. Let us distin- 
guish four of them: 

• preheating, with the duration bt\ ^ lOO/m, 

• a short violent stage at the end of preheating when 
non-linearity and non-perturbative effects are dom- 
inant; chaos is onset, erasing the details of the ini- 
tial conditions, on a timescale bti ~ lO/m, 

• the stage of turbulence of classical fields, with the 
saturated spectrum cascading both towards UV 
and IR modes with duration btj,^ longer than the 
two previous stages, 

• the last stage of proper QFT thermalization, with 
particle fusion/offshell processes; it is characterized 
by the timescale (5^4, much longer than the previous 



Traditionally the beginning of the thermal history is 
related to the moment (5^4, which may be very long; as 
a result of redshift, the reheating temperature can be 
rather low. 

Let us reconsider this attitude. While the question 
of when the ultimate thermal equilibrium will be estab- 
lished and to which value of it corresponds is still very 
interesting, we will argue that important physics which 
constrains the model of inflation and interactions takes 
place long before full thermal equilibrium is completed. 

Indeed, Figures EHOl show that the modes at physically 
interesting scales, up to fc '--^ tens of m ~ 10^^ GeV, are 
quickly excited towards a saturated distribution. This 
raises several subtle and interesting questions. Usually, 
the relativistic particles embedded in the thermal bath 
are brought into thermal equilibrium after the relaxation 
time 1/crn, where a is the cross section of the processes 
relevant for thermalization, and n = J (Pkuk is the num- 
ber density of particles. This estimate works in the di- 
luted gas approximation. Indeed, consider the vertex 
0X^1 leading to inflaton decay into quanta of x- The 
particle number density enters in the coUisional integral 
in the combination 



'pi 



(1 + n^2) (1 + n^,) - <3 (l + <i) (23) 



(pi is the momentum of the inflaton quantum, while 
P2.3 are the momenta of the two quanta of x entering 



C. Physics before thermalization 

One of the most important parameters in physics of 
the early universe is the highest temperature of the hot 



^ It is well known that if you assume a slow, perturbative inflaton 
decay and a quick thermalization of the decay products, a sub- 
dominant thermal bath will be present while the inflaton is still 
decaying. That will have a temperature muich higher than Tr, 
defined at the moment when the thermal bath starts dominating. 
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in the process). The above estimate for the thermahza- 
tion timescale assumes that the inverse process XX ~* 4' 
(described by the second term in (j2SJ) is irrelevant, and 
that all the occupation numbers are much smaller than 
one, so that 1 + rip ~ 1 . In particular, this last assump- 
tion is not valid until the very end of the thcrmalization 
stage. 

Therefore, for large Uk we expect that particles are 
dragged into effective thermal equilibrium faster, due to 
the stimulated interaction. We have observed this fact 
also in previous numerical lattice simulations ,12,]. 

Another relevant effect, already discussed above, is 
that large fluctuations of fields can contribute to the ef- 
fective masses of particles, which can become effectively 
heavy. This leads to the blocking of some processes which 
may lead to a faster thcrmalization, until fields fluctua- 
tions are diluted by the expansion. 

Finally, a very relevant question is the one of parti- 
cle production from the nonthermal but highly excited 
distributions. In particular, one should verify that dan- 
gerous relics are not overproduced at this stage. To see 
this, let us discuss gravitino production from the decay 
products of the inflaton. In models where supersymme- 
try is broken gravitationally, only a very small number 
of gravitinos can be tolerated, 1^/2 = n-ii^js ^ 10~^^ 
(the exact value being dependent on the gravitino mass) 
|19|^. The typical approach is to compute the produc- 
tion of gravitinos only after the thermal bath has formed. 
The dominant processes are 2 — > 2 scatterings with only 
one gravitino as outgoing particle (and, hence, only one 
gravitationally suppressed vertex, while the other vertex 
is typically a gauge interaction). The rate of gravitino 
production follows from the Boltzman equation 

^^ + 3HN-,^2 = {'j\v\)NtNt , (24) 
at 

where Nt ~ is the abundance of the (MSSM) degrees 
of freedom in the thermal bath, and where the average 
cross section is of the order {a\v\) ~ a/M^ , where a is a 
gauge structure constant. The right hand side decreases 
very rapidly with time, and the production at the highest 
possible temperature (that is, at Tr) dominates. Hence, 
only the production in the first Hubble time is relevant, 
and the final result can be estimated as 

^3/2 - ^ X H^^' , (25) 

where Hr ~ lOT^/Mp is the Hubble parameter at T = 
Tr ■ This leads to the gravitino abundance 

I3/2 - 10-^ ^ , (26) 



^ A more stringent limit app lies if the decay has a significant 
hadronic branching ratio l26l . 



and to the well known bound Tr ^ 10 GeV on the re- 
heating temperature. 

This standard computation of gravitinos generated 
from the thermal bath neglects all the production which 
may have taken place during the thcrmalization stage. 
The underlying idea is that, during reheating, most of the 
energy is in the inflaton condensate, where the quanta 
(loosely speaking) do not have momentum, and hence 
cannot scatter to produce gravitinos. While this may be 
true for a slow (perturbative) inflaton decay, this is cer- 
tainly not the case for preheating/rescattering, where the 
distributions of the decay products form right after the 
end of inflation. These distibutions may be responsible 
for significant gravitino production, even before complete 
thcrmalization has taken place [27|. 

The computation is now more difficult than in the ther- 
mal case 19]. First of all, unlike the thermal computa- 
tion, the result is expected to be model dependent. Sec- 
ond, turbulent nonlinear processes are dominant right 
after preheating, so that it is possible that processes in- 
cluding more vertices give a sizeable contribution to the 
amount of gravitinos produced. However, we can obtain 
a conservative estimate by considering only 2^2 tree 
level processes also in this case. For definiteness, let us 
take the model ([[J, Q, with 5^ = 2.5 x 10"'^ (the value 
that we have studied above), and assume that the x field 
produced at rescattering has a trilinear vertex with its 
femionic superpartner x and another fermion z. For in- 
stance, this could be a gaugino, if x has gauge interac- 
tions. However, let us consider here a generic vertex with 
coupling constant h and the interaction xX ~^ 4^3/2 
mediated by the fermionic partner of x- (As a technical 
point, we consider production of transversal gravitinos, 
since the nature of the longitudinal component changes 
with time). 

We are interested in evaluating the production at the 
time ~ 120/rn, when the (physical) number density 
of quanta of x is maximal. We first discuss whether the 
process we are considering is kinematically allowed. At 
the time , the non-gravitational vertex provides an ef- 
fective mass for fermions of the order (see Fig. I1U|) . 

(U) - h V(x2 (U)) - 5 • 10^ /im (27) 

The distributions of x obtained from the lattice have a 
typical comoving momcnum of the order of 10 m , which, 
at the time , corresponds to the physical momentum 
p* ~ 0.5 m . Hence, ii h < 10~^ , most of the quanta of 
X will be able to produce gravitinos through the process 
we are considering. For higher values (for example, if the 
second vertex is a gauge interaction), only the few quanta 
of X with very high momenta can contribute to the pro- 
duction, so that this process is signifi cantly suppressed. 
However, as the universe expands, vO?) decreases as 
a~'^/^ , while the physical momenta p cx decrease 
less. Therefore, more and more quanta will have a suffi- 
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ciently high momentum to produce a tp3/2~z pair.^ For 
simphcity, h <, 10"'^ is assumed in the following discus- 
sion; we do not expect that the following conclusions will 
significantly change if a stronger coupling is considered. 

By proceeding as in the thermal case, we can estimate 
the number of gravitinos produced in the first Hubble 
time after . (Again, the numerical values are taken 
from the lattice results; (i*)^ denotes the physical 
number density of the x quanta at the time t, .) 

N3/2 (U) -^N^ {t*f H iU)-' ^ 10'' hm^ . (28) 

p 

For sufficiently high h , this value is actually incorrect, 
since it exceeds Pauli blocking. Assuming that the whole 
Pauli sphere up to physical momentum is filled gives 

Ny^iQ-O.lm^ (29) 

The number density of produced gravitinos is then the 
smaller between eqs. (I28II and H29(l . according to the value 
of ft. . In the following, we assume h > 10^* , so that Pauli 
blocking is reached and the number density of gravitinos 
is given by eq. f^ . 

The gravitino abundance is obtained by dividing the 
gravitino number density by the entropy. At the time 

, the distributions of the decay products are far from 
thermal (being concentrated in the IR), so that their en- 
tropy is smaller than that of a thermal bath with the 
same energy. The entropy density can be computed from 
the lattice simulations to be approximately s — 10'*™^ 
at this stage |2^. However, the (comoving) entropy is ex- 
pected to increase due to thermalization. If we convert 
the actual particle distribution instanteneously into the 
thermal distribution, the entropy would reach the max- 
imal value Smax ~ 10^ allowed by energy conserva- 
tion. The difference between the actual entropy density 
and Smax gives rise to the dilution of the gravitino abun- 
dance during thermalization. (Further sources of dilution 
require extensions of the minimal model.) 

Therefore, the conservative estimate for the gravitino 
abundance is obtained by dividing the gravitino number 
density by 

p(t,)'/*- lO^m^ . (30) 

where p is the (physical) background energy density at 
the time . Hence, the ratio 

- - 10-^ (31) 



* This scaling argument does not include the effect of the thermal- 
ization, which is also increasing the p/mg ratio. Indeed, since 
thermalization proceeds through particle fusion it has the effect 
of both increasing the typical momenta of the distributions and 
of decreasing the effective masses. 



is about 6 orders of magnitude higher than the allowed 
limit. 

This is a serious concern, which has to be addressed 
in the model building of inflation. It can be resolved by 
introduction of extra radical assumptions. For instance, 
this gravitino abundance can be decreased by an entropy 
dilution. This happens if a massive species ip is produced 
at reheating, and it dominates the energy density of the 
universe for some time. When it decays, it leads to a 
significant amount of entropy, which dilutes the parti- 
cle pieces produced directly at the inflaton decay. Con- 
cretely, the entropy dilution leads to the decrease of the 
gravitino abundance ~ ^■ijil f?^^ ■ Supposing that 
the universe is matter dominated between the times t\ 
and t2 , the ratio N^/2l P^^"^ is decreasing as (01/02)'^^'* 
(fli being the scale factor at the time ti). This entropy 
dilution could also occur through a secondary stage of 
inflation "29^. 

In the specific toy model we have considered, a decrease 
of six orders of magnitude does not seem to be realistic, 
so that the production of gravitinos in this model ap- 
pears to be a serious concern. However, it is important 
to remember that the result H31|l is very sensitive to the 
specific model considered, and that a case by case calcu- 
lation is necessary. 

D. Modulated Fluctuations from Preheating 

One of the motivations for studying the EOS after pre- 
heating is related to the theory of modulated cosmolog- 
ical fluctuations. We briefly review the idea as it was 
presented in 0| , and discuss the implications of our 
findings on the evolution of the EOS for this mechanism. 
We leave the explicit calculations for subsequent publi- 
cation. 

Suppose that the coupling depends on some modu- 
lus field z, (as is typical in string theory) 

9^^9\z), 6g^^^5z. (32) 
az 

During inflation the modulus field z can be light, with a 
mass smaller than the Hubble parameter during inflation 
H ^ 10^^ Gev. Then, large scale fluctuations 5z with an 
almost scale free spectrum are inevitably generated from 
inflation. The wavelength of fluctuations exceeds the size 
of the causal patch H~^. As we have discussed above, the 
EOS varies very quickly already at the beginning of the 
thermalization stage. The variation is a non monotonic 
function of the coupling g , see for instance Figs.n]and|21 
Due to the large scale spatial fluctuations of the coupling 
Sg^, the change in the EOS proceeds with a slight time 
shift in different Hubble patches. This results in the gen- 
eration of (almost) scale free scalar metric cosmological 
perturbations. 

To perform accurate calculations of the amplitude of 
modulated fluctuations, one needs to know the exact 
character of the transition from inflaton fleld to radiation. 
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In the simplified picture of perturbative refieating, when 
the inflaton energy decays exponentially as ~ e"'"*, 
calculations of modulated fluctuations |3| are based on 
linearization of the inflaton decay rate F, which is lin- 
early dependent on the coupling . However, as we have 
remarked, the decay of the inflaton field is typically more 
complicated, and the transition towards RD occurs very 
quickly (well before thermalization completes). In addi- 
tion, we stress that the time of the transition is a non 
monotonic function of . Finally, although we focused 
only on the direct coupling of the inflaton to matter, 
the details of thermalization (and of the evolution of the 
EOS) depend on the other interactions among the decay 
products; the strength of these interactions are likely to 
be dependent on some modulus field as well, and they 
can also modulate cosmological perturbations. All these 
features have to be taken into account for precision com- 
putations of modulated fiuctuations. 



V. SUMMARY 

In this paper we studied the out-of-equilibrium nonlin- 
ear dynamics of the fields during and after preheating. 
This continues earlier works in [3) llsL IT^ . We use lat- 
tice numerical simulations of fully non-linear dynamics 
of interacting classical fields. In this investigation, we 
computed the evolution of the equation of state during 
preheating and the early thermalization stage. Immedi- 
ately after inflation, the EOS is the one corresponding 
to inflaton domination. For a massive inflaton, we have 
matter domination, characterized by p = 0. After a very 
short time, the inflaton transfers its energy into inhomo- 
geneous modes. At the first stage (preheating) inhomo- 
geneous field fluctuations are copiously produced in the 
regime of parametric resonance. Then the modes rescat- 
ter to redistribute the energy by cascading in the phase 
space. This process of classical wave turbulence can be 
rather long. This leads to a final stage, where the occupa- 
tion numbers are small and the classical approximation 
breaks down. Thermalization finishes when the quantum 
fields reach the ultimate thermal equilibrium. If all the 
participating fields are light, the final thermalized state 
corresponds to an EOS of radiation, with p = ie. 

The transition between the matter and radiation dom- 
inated EOS is hence typically supposed to occur on a 
(relatively) very long timescale. What we found in our 
numerical simulations is very different from this naive ex- 
pectation. In fact, the EOS jumps from matter domina- 
tion to (almost) radiation domination immediately after 
the first stage of preheating, i.e. in a couple dozen in- 
flaton oscillations (10~^^ sec each), long before thermal 
equilibrium has been established. In other words, the 
macroscopic EOS is close to the one of radiation domi- 
nation, while the microscopic state is far from thermal. 
This result is similar to the conclusion of ^3 regarding 
0{N) theory in flat spacetime. 

The sharp change of the EOS we observed can be very 



conveniently related to the moment of energy transfer 
from the inflaton to the inhomogeneous radiation. As we 
mentioned, it is a very fast process, as is well-known from 
the theory of preheating, which takes place in only about 
10"'^^ sec. Moreover, the time at which the transition oc- 
curs depends only weakly on the coupling g^ between the 
inflaton and the other flelds. The exact dependance is 
non-monotonic, and the transition time oscillates around 
^ 10~^^ sec (as set by preheating) for a very wide range 
of couplings. This result is drastically different from per- 
turbative inflaton decay, where the transition timing is 
proportional to the coupling. 

This very quick evolution has several consequences, 
and leads to a number of issues which deserve further 
investigation. We stress four of them. 

i) The first is related to the generation of modulated 
cosmological fluctuations from preheating. It is conceiv- 
able 0, 13, 01 that the strength of the inflaton interac- 
tions with its decay products is a function of some moduli 
flelds as expected in string theory. Fluctuations of the 
light moduli (with mass smaller than 10^'^ GeV) are in- 
evitably generated during inflation. This results in spa- 
tial variations of the couplings at large scales, well beyond 
the size of the causal (Hubble) patch at the end of infla- 
tion. Hence, the transition from inflation to radiation 
domination occurs at slightly diffferent times in different 
Hubble patches, giving rise to scalar metric fluctuations. 

To calculate metric fluctuations from modulated per- 
turbations, one needs to know the exact evolution of the 
EOS, from the matter domination (inflaton) to the radi- 
ation domination (light decay products) stage. We found 
that the transition of the EOS occurs sharply in a step- 
like manner. Therefore one can introduce a time-like (but 
spatially rippled) hypersurface which divides two cosmo- 
logical regimes, with p = and p — ^e. Metric fluctu- 
ations can be calculated using GR matching conditions 
across this hypersurface, similar to how it was done to 
calculate modulated fluctuations in hybrid inflation |1C| . 

ii) The second application of the EOS history is related 
to the formula which links the number of efoldings during 
inflation N and the present day physical wavelengths of 
fluctuations in terms of the logs of their momenta logfc. 
Indeed, the cosmological evolution of the scale factor a (t) 
is ultimately defined by the EOS. Refining the A'' — log k 
relation is mandatory to put precise observational con- 
straints on infiationary models. In this paper we notice 
that a (practically) immediate transition from inflation 
to radiation domination is a common feature of inflation- 
ary models with preheating, either parametric resonance 
or tachyonic ones. Equation H21|) . with our result H22|l . 
refines the A^ — log k relation. 

iii) The next note is related to the model building of in- 
fiaton interactions. In the majority of papers discussing 
infiaton decay, interactions with bosons are considered in 
the form of four- legs vertices of the type g'^ffP'x^ . These 
interactions are the dominant ones at the early stages of 
preheating, when the amplitude of the infiaton is large. 
(We do not consider here nonrenormalizable interactions. 
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which - if present - could be more important.) However, 
they do not lead to a complete inflaton decay. A com- 
plete decay requires vertices with only one inflaton, as 
the Yukawa interaction (jjipip wih fermions, or the three- 
legs bosonic interaction ct^^x^ (where cr is a mass scale). 
As we have argued, these trilinear interactions are very 
natural in supersymmetric models. Supersymmetry has 
the additional advantage of protecting the flatness of the 
inflaton potential from too large radiative corrections, 
which arise from the coupling of the inflaton to the light 
degrees of freedom. 

iv) Finally, the early evolution of the light degrees of 
freedom from the decay of the inflaton can have very 
important consequences for the generation of particles 
far from thermal equilibrium. A particularly relevant 
application is the possible (over) production of graviti- 
nos from the decay products of the inflaton. Although 
this is a perturbative generation, it can be very effective, 
since it occurs at very early times. Traditionally, the 
limit Tr ^ 10^ GcV is due to gravitinos produced only 
after the x quanta have thermalized. However, a non- 
thermal distribution can also be responsible for gravitino 
overproduction . These inflaton decay products arise 
very rapidly, with an energy density much higher than 



(lO^ GeV) ; although their initial distribution is far from 
thermal, we have seen that they can be responsible for 
a significant gravitino production, which in some cases 
exceeds by several orders magnitude the allowed bound. 
Clearly, this is a less universal bound than the standard 
one, since it depends on the details of how the thermal 
bath is produced. (For instance, the gravitino abundance 
can be diluted by entropy release during thermalization.) 
In addition, the gravitino problem is avoided altogether 
if the gravitino mass is significantly smaller (as in gauge 
mediated supersymmetry breaking) or larger (as in the 
simplest versions of the KKLT model [3^ ) than the elec- 
troweak scale. Still, this possible production must be 
addressed in many inflationary models. 
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